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Abstract: On account of its application to the present and future analysis of the 
virtual photon correction to the hadronic properties by means of lattice QCD sim- 
ulation, we initiate the study of the finite size scaling effect on the QED correction 
using low energy effective theory of QCD with QED. For this purpose, we begin 



with formulating a new QED on the space with finite volume. By adapting this 
formalism to the partially quenched QCD with electromagnetism, we explore the 
qualitative features of the finite size scaling effect on the electromagnetic correction 
to the masses of pseudo-Goldstone bosons. 
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1. Introduction 

The recent progress in the lattice QCD simulation enables us to perform "measure- 
ments" for hadronic properties through the world with more and more realistic QCD 
realized in the computer. It may not be a far future that the lattice QCD simulation 
becomes a method by which precise measurement of nonperturbative dynamics of 
QCD is reached. 

One of the important applications of the lattice QCD simulation would be the 
determination of quark masses. The light quark masses have been determined using 
dynamical quarks with two flavors [|, g |, § § § and (2 + 1) flavors 0, § |, |nj. 
The forthcoming precise measurement through the lattice QCD simulation reminds 
us that quarks are electrically charged. All of the hadronic properties thus suffer 
from electromagnetic (EM) radiative corrections. Because the EM interaction is 
other source of explicit breaking of isospin symmetry than the difference between 
the masses m u , m d of up and down-quarks, the determination of m u — m d requires 
us to grasp the size of EM correction quantitatively at the hadronic level. Among 
the references [0 |, g |, |, % % g |, 0, Ref. f§ is the only work that presented 
the values of m u , m^ each, but it does not seem to make the data fit with a function 
parametrizing the EM splitting in the kaon masses as an undetermined constant. 
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The attempt to incorporate EM correction to the pseudoscalar meson masses in 



the context of lattice simulation was first done by Duncan et g|. ||ll|| by simulating 
QED put on the lattice. By adopting the same technique to incorporate QED cor- 
rection, the QED correction to the meson masses has been investigated with use of 
dynamical domain wall fermions with two-flavors [12| , and in the quenched approxi- 
mation with the renormalization group improved gauge action [|T^]. Ref. (JTJJ employs 



another method to calculate the leading-order EM correction to Am 2 = rn 2 + — m 2 

O 7T 7T ' 7T U 

in the two-flavor overlap fermion simulation. This work relies on the formula Jl5 
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where D^ u (k) = rf" jik 2 + ie) is the photon propagator, /„. (~ 92 MeV) is the 
pion decay constant, VT = qT "^^ and A a = qT a ^f iy ^f 5 q with the normalization of 
SU(2) generators tr (T a T b ) = \5 ah . (^0)q CD denotes the expectation value of the 
operator O with respect to QCD. The lattice simulation calculates the two correlation 
functions appearing on the right-hand side of Eq. ( |1 . 1| ) and attempts to get the pion 
mass difference in the chiral limit [ |T4|j . 

Though our interest is the nonperturbative dynamics of QCD on space with 
infinite volume, the simulation has to be done in the virtual world with finite volume. 
As quarks with colors and electric charges must live in such a finite volume, it 
is inevitable that QED correction suffers from finite size scaling effect no matter 
what computational method one may choose. It is plausible that two pseudoscalar 
mesons in a common isospin multiplet share the same QCD finite size corrections 
[0, 0, 0, [Zg, gg, £2} [2§ £§ £Z§ [ZB|, 0, 0. In addition, electromagnetic force 



is long-ranged. Therefore the finite size scaling effect regarding QED will dominate 
the finite size corrections in the EM splittings. Unless it is adequately quantified, 
ignorance on the QED finite size scaling effect becomes another source of systematic 
uncertainty in m u — m^ derived using lattice simulation. 

Thus far, there are diametrical opposite views on the relevance of QED finite size 
scaling. The authors in Ref. [|IB[] performed the direct (quenched) lattice QCD mea- 
surement of the EM splitting in pion masses using two different sizes of four- volumes, 
T x L 3 = 24a x (12a) 3 and 24a x (16a) 3 , where a is the lattice spacing. They observed 
no significant difference for the EM splitting measured in two volumes and concluded 
that a linear size L = 2.4 fm is sufficient to compute the EM splitting without cor- 
recting the measured values within available statistical uncertainty. Contrastingly, 
Refs. |H| , |1~2]| estimated the QED finite size scaling according to the one-pole satu- 



ration approximation [IS, ^ ^] to Eq. ( |1 . 1| ) with the momentum integral replaced 
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with the sum 



a 2 1 3a 2 1 1 v^ m\rm 2 A 

fee(f 4 -{0}) viy A> 



where my ~ 770 MeV, m^ ~ 970 MeV is the mass of A\ in the chiral limit, and 

2n 2n 

k = (h , h, k 2 , h) k e — z , kj e — 



koe—z^je—z}, (1.3) 



is a lattice on the Euclidean four-momentum space. For the lattice geometry TxL 3 
32a x (16a) 3 with 1/a ~ 1.66 GeV in Ref. |2j, Eq. Q leads 



Aml\ VMn (T, L) 



Am ^L,VMD(°°' °° 



which is not negligible. 

Under such circumstances, we investigate the finite size scaling effect on the 
QED contribution to light pseudoscalar meson masses using the chiral perturbation 



theory including electromagnetism J311 R3, 33, 34, 35]. Though the finial quantitative 



determination of the finite size correction must resort to the first principle calculation 
as in Ref. ]T3f , any qualitative understanding as shown here will help one to make 



extrapolation to infinite volume in the future simulation study. 

This paper actually consists of two parts. In the first part, Sec. |2|, we give a 
formulation of QED on the space with finite volume. This task is necessary because 
a single classical charged particle does not satisfy the equation of motion, say, the 
Gauss' law constraint in the finite volume QED obtained via the ordinal compacti- 
fication procedure. This is intuitively understandable as the electric flux emanating 
from a nonzero charge finds nowhere else to go. We thus begin with defining such 
a QED that accommodates a single charged particle on the space RxT 3 , where T 3 
is three-dimensional torus corresponding to the compact space. The second part, 
Sec. |3|, utilizes this QED to study the finite size scaling effect on the electromag- 
netic splittings in pseudoscalar meson masses. For this purpose, the calculation is 
performed in the framework of the partially quenched chiral perturbation theory 
[|36| , [37], |38], |39j including the electromagnetism |33| in view of its practical applica- 



tion to the actual lattice simulation . Sec. f| is denoted to discussion and conclusion. 
Appendix [A] collects the formulae for the basic sums which appear in the evaluation 
of finite size scaling effect. 

2. QED in finite volume 

The aim of this section is to present a new QED on the space with finite volume, 
which allows us to investigate the properties of a single charged particle. We first 
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clarify in Sec. |2.1| the problem itself that confronts us in the QED obtained by the 



ordinal compactification procedure. In Sec. [2.2| , we define a new QED in finite volume 
and explain how it solves this problem. Throughout this paper, the topology of the 
space is the three-dimensional torus T 3 = S 1 x S 1 x S 1 with a common circumference 
L for all S 1 . A point on T 3 is thus specified by the coordinates x = (x 1 , x 2 , x 3 ) 
obeying periodicity x J = x J ' + L (j = 1, 2, 3). As in the analysis of finite size scaling 
of QCD [jnj |T7], [TS|], the temporal direction t = x° is taken to be infinite, tGK, for 
single particle states to develop poles in the energy space. We adopt the convention 
r)fj, u = diag (1, — 1, —1 — 1) for the signature of the metric. 

2.1 problem 

The most familiar procedure to construct the corresponding theory on T 3 is to impose 
periodic boundary conditions on all fields. The electromagnetic theory obtained via 
such a naive Z 3 -orbifolding procedure 1 is referred to as QEDjz here. In QED Z 3 
the gauge potential A^x) [x v = (t, x)), in particular, obeys periodic boundary 
condition in every spatial direction. This form of boundary condition is motivated 
for the practical reason that the periodic or anti-periodic boundary condition is 
imposed along spatial directions for the available lattice QCD configurations and the 
conservation of quantized momenta at every QED vertex requires that every spatial 
component of photon momenta be an integer multiple of ^. The action for the gauge 
kinetic term in QED^s takes the usual form 

S y = J dt J ^ d 3 x f-i Fp, fA , (2.2) 

with the field strength F^ v = d^A u — d v A^. 

To elucidate the problem from the practical point of view in our context, we go 
back to Eq. (|1.2j ) for the EM splitting in the pion mass squared in the effective theory 
including the vector and axial-vector mesons. The application of the Z 3 -orbifolding 
procedure to the effective Lagrangian including the vector and axial-vector mesons 
in Ref. |30| immediately leads the same form of the Lagrangian written in terms of 
the spatially periodic fields. This effective theory gives the expression for the EM 
splitting in Eq. (|1.2|) with T — ► oo 

A -.L,vmd (°°> L) = - (4.) J — - ^ fc2(fc2 + m2)(fc2 + mi) , (2-3) 

ker 3 



1 The word naive means here that the Fourier modes A^t, k = 0) in Eq. ( |2.8| ) are treated as 
basic variables. We recall that the Wilson lines is appropriate variables. The correct Z 3 -orbifolding 
would involve the integration over Wilson line Uo(t) ~ exp ie J dt j^ Ao(t, 0) which leads to the 
constraint that the charges on every three-dimensional hypersurface should vanish in total 

-^|d 3 x/(i,x) = 0. (2.1) 

We would like to see what happens in the naive procedure. 
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where 

= (k\k 2 ,k 3 ) We— z| . (2.4) 

We can see that the quantity ( |2.3| ) has infrared (IR) divergence coming from the 
contribution of k = 0. The expression ( |2.3| ) can be interpreted as the sum of the 
contribution of Kaluza-Klein modes in one- dimensional field theory. The IR diver- 



gence appearing in Eq. Q2.4Q is attributed to that of a massless mode from this 



point of view. It should be reminded that Am^\ a VMD in infinite volume is IR-finite. 
Therefore, lim Am1\ .„„, (oo, L) ^ Am^l ,„„,. The theory deducing Eq. (12.3) 

L— >oo " I a, VMD v 7 ' "la,VMD ^ o ± m \j 

is one of QED^s with the content of matter fields and the action concretely specified. 
The same pathology always emerges in every QED Z 3 theory irrespective of the details 
of the matter fields and the action. Thus QED Z 3 cannot be used to study the finite 
size scaling effect on the EM splittings. 

The origin of the above pathology will be traced back to the inconsistency of 
a single charged particle with the classical equation of motion. Though it may be 
well-known, the detail observation of this point will help us to grasp the essence of 
our new QED. Due to the periodicity of the gauge potential, the electromagnetic 
current j M (x) = SSmaxter/SA^x) derived from the matter part Scatter of the action 
is also periodic along the spatial directions. The current j fJ '(x) is assumed to be 
conserved under the equations of motion derived from the variation of matter fields. 
The classical equation of motion derived from the variation of A^x) is hence 

d v F* u '{x)=f(x). (2.5) 

This contains the Gauss' law constraint 

V • E(x) = p(x) ■ (2.6) 

where the electric field E^(x) and the charge density p(x) are given by E^(x) = F°i(x) 
and p(x) = j°(x), respectively. For simplicity, we consider an infinitely heavy charged 
particle. When it is at rest initially, the charge density and the profile of electric 
field are both constant in time 

V-E(x) = e<5 3 (x). (2.7) 

The inconsistency appears when both sides of this equation are integrated over the 
whole T 3 ; the left-hand side vanishes while the right-hand side does not. Likewise, 
we can see that any single charged particle cannot live on T 3 . 

2.2 new QED onlxT 3 



The observation in Sec. 2.1 shows that the IR divergence in Eq. (|2.3|) is a mani- 



festation of inconsistency of a single charged particle with the classical equation of 
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motion. The aim of this section is to introduce an alternative QED in finite volume 
that solves this problem for the sake of our study of the finite size scaling effect on 
the EM splitting. 

In the space T 3 , three-momenta take discrete values in T 3 in Eq. ( |2.4|) . Accord- 
ingly, the gauge potential is decomposed in the Fourier series with respect to the 
spatial dimension 



Mt> 



i £ e *-M„(t, k) . 



:2.8i 



ker 3 



The new QED on K x T 3 , referred to as QED^, is the theory without variables 
Afj,(t, k = 0) ab initio. In other words, we do not incorporate the Wilson lines 
U,j,{t) ~ exp ie f dt j^A^(t, 0) as dynamical variables. In terms of such a gauge 
potential, the action of pure electromagnetism is given by 5 7 in Eq. (|2.2| ). Below, we 
observe the various features possessed by QED^. 



First, we see how QED^ solves the problems in Sec. 2J . The equation of motion 
is modified as follows. Since the modes A^it, 0) are absent, the variation of the full 
action S = S 7 + Scatter with respect to the gauge potential becomes 



SaS 



dt 



-oo 

oo 



cfx 



.-SF^ + SA.f 



dt J2 l s M^ k ) 

kef; 



where 



x / ^xe lk - x (-^F^(t, x) + fit, x)) 



r' 3 = r 3 - {0} 



[2.9) 



(2.10) 



As the result, the extremum condition gives the equations of motion only for k ^ 0; 
/"t/ 3 xcos(k-x){^F^(t, x)-f(t, x)} = (k^0). (2.11) 



The inconsistency seen in Sec. |2J] is therefore circumvented. We also note that 
Eq. (|2. 11|) instead of Eq. (|2.6| ) no longer yields the equality between the charges 
contained in a domain V in T 3 and the electric flux penetrating the surface dV. 

The quantum field theory will be defined in the path integral framework by 
defining the measure of the gauge potential as usual. Let SA^t, k) be an infinitesimal 
variation of the mode A^t, k) (t 6 M, ke T' 3 ). The functional measure is defined 
corresponding to the norm in the space of gauge configurations in QED^ 



\\SA\f 



f dtJ2 5A ^ k)8A»{t, k) 
kef; 



(2.12) 
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As we will see later this norm will turn out to be gauge- invariant. From the point 
of view of field theory on R, A(t, 0) are massless fields. Due to the absence of these 
modes in QED^, the IR divergence as in Eq. Q2.4|) no longer appears. 

We next observe the gauge symmetry of QED^. The transformation of the gauge 
potential is given by 

A lt (x)»A' li (x)=A tl (x)+d ll A, (2.13) 

while a matter field $(x) with charge Q$e is transformed as 

$0) ^ $'(z) = exp[iQ<s, e A(x)] $(x) . (2.14) 

Here we assume that there is at least one matter field which has the minimum charge 
e. The easiest way to identify the full gauge symmetry is to look at the full gauge 
symmetry of QED Z 3 first. There, the general form of A(x) that keeps the periodic 
boundary condition for the matter fields takes the form 

2tt 3 
A(x) = A P (x) + — ^mjx* (QED z a) , (2.15) 

i=i 

where rrij G Z, and Ap(x) is periodic along T 3 . 

Denoting the Fourier components of Ap(t, x) by Ap(t, k), the gauge transfor- 
mation for the gauge potential becomes in the three-momentum space 

A'fi, k) = Aj{t, k) + (ik j A P {t, k) + L 3 5 k ,o ^ mA (QED z3 ) . (2.16) 

In our new QED, A(t, 0) no longer exists. It thus turns out that there are no 
redundancy corresponding to m G Z 3 and time- dependent but spatially homogeneous 
part of Ap(x), that is, 

4(£, k) = A s (t, k) + % k j A P (t, k) (QEDJ , (2.17) 

where 

d t A P (t,0)=0. (2.18) 

One can readily see that the set of all functions that fulfill Eq. ( |2.18| ) forms an abelian 
group. From Eq. (|2.17| ), it is also easy to see that this gauge group is exactly the 



redundancy that allows us to take the Coulomb gauge fixing condition 

djAjit, x) = 0. (2.19) 

We recall that in QED^3 a residual gauge symmetry survives even after imposing the 
condition fl2.19| ), which should be fixed by the additional condition A (t, 0) = [|T^]. 



-7- 



In contrast, the condition Q2.19 ) suffices to fix redundancy in QED^ leaving only the 



global symmetry. Obviously, the norm (|2.12|) in the space of gauge configurations is 



gauge-invariant. We can also make a BRST complex by introducing the ghost fields 
corresponding to the gauge parameters of the form ( |2.18| ) in the standard manner 

[0- ' 

We close this section with a few remarks. We consider the scattering of a charged 

particle and its anti-particle in the center of mass frame, leaving aside the issue 

whether the scattering may not be well-defined in the presence of a long-ranged 

force in finite volume. Due to the absence of the modes A^(t, 0), the s-channel 

process mediated by a single virtual photon does not occur. However, we recall that 

Lorentz invariance, in particular, the symmetry related to the Lorentz boosts, is 

explicitly violated on the space K x T 3 . Thus, once we consider the collision say, of 

an incident charged particle with three-momentum (p+ ^-)e x , where p G ^Z and e^ 

is a unit three- vector along x-direction, and an anti-particle with momentum (— p)e x , 

the s-channel process occurs. In the limit L — » oo, the cross section will approach to 

that in the center of mass frame in infinite volume. 

Secondly, since the equation of motion ([2.11|) is not written locally, QED^ seems 



to possess somewhat non-locality In fact, it is not still clear at the present stage 
whether this is actually the case, and then whether another pathology appears in 
QED^. As will be demonstrated explicitly in the subsequent section, even if non- 
locality is present, it is so mild that the structure of ultra-violet (UV) divergence 
in QED^ remains completely the same as in QED in infinite volume. This feature 
is contrasted the situation in noncommutative field theory Hi], 02], H3]; the the UV 



structure of a noncommutative field theory differs significantly from that of the com- 
mutative counterpart due to the hard non-locality 

3. Finite size scaling in meson mass 

Now we apply QED^ to the study of finite size effect on the EM correction to the pseu- 
doscalar meson masses in the chiral perturbation theory including the electromag- 
netism [fTi], |32], |33|, Q. Looking at the practical application to the lattice simulation, 
we adopt the partially quenched chiral perturbation theory [J3(| [37], |38], [3^| including 



the electromagnetism J33| and compute the leading-order finite size correction to the 
EM splitting in this theory. For that purpose, we begin with summarizing our nota- 
tions for partially quenched chiral perturbation theory including electromagnetism to 
the next-leading order. We derive the formulas for the next-to-leading order correc- 
tion to the pseudoscalar meson mass in finite volume, and evaluate them numerically 
to investigate the finite size correction to the EM splitting. 

3.1 partially quenched chiral perturbation theory with electromagnetism 

The next-to-leading order corrections to the off-diagonal meson masses has already 



been computed in Ref. ]33|]. However, the expression written by the momentum in- 
tegrals are needed in practice for the study of finite size scaling effect. To derive 
such an expression in Sec. |3.2| , we fix the notations, in particular, of the low-energy 
constants at the next-to-leading order for the subsequent calculation, and dictate 
the free meson propagators necessary for the one-loop calculation. To take the ap- 
plication to the partially quenched system with two-flavors, the super-trace of the 
EM charge matrix Q is not assumed here to vanish, unlike in Ref. |53| . There then 
appear more local terms at the next-to-leading order (i,e., 0(p 4 ), 0(e 2 p 2 ) and 0(e 4 )) 
than those found in Ref. |33 |. We list up all of them by drawing upon Appendix of 



Ref. [32|| where one-loop UV divergences were computed for generic flavor number 
Np in the unquenched chiral perturbation theory including electromagnetism. After 
that, we compute the UV divergences that should be absorbed by the coefficients of 
these local terms. 

In what follows, all the fields and parameters are written in the "flavor" basis 

O. = (qY, ■■■ , qN v , gf, -•• , $ g , gi, ••• , 9n v ) , (3.1) 

where qf (r = 1, ■ ■ • , N$) denote the sea quark fields, q^ (a = 1, • • • , Ny) the 
valence quark fields, and g a (a — 1, ■ • • , Ny) the ghost quark fields f36fl . The chiral 
symmetry is a graded Lie group G = SU(N S + N V \N V ) L x SU(N S + N V \N V ) R . It 
breaks down spontaneously to its vector-like subgroup H = SU (N s + N v \N v )v ■ The 
associated Nambu-Goldstone bosons are represented by an (Ns + Nv\Ny) x (N$ + 
Nv\Nv) supermatrix II. Using 

u[U(x)] = exp (i iM-) , (3.2) 

II transforms nonlinearly under (gi, g R ) G G through 

u[U] ^ u[W] = g R u[U] h((g L , g R ); n)t = h((g L , g R ); n)u[Ii]g{. (3.3) 

where h((gL, g R ); II) e H. We follow the convention for the chiral Lagrangian which 
can be read off from Refs. |33|, |44| with minor modification. 

The external fields i? M (x), L^x) that couple to the right-handed and left-handed 
chiral components of vector currents are incorporated in the partially quenched QCD 
action for the purpose of calculating the connected Green functions of vector and 
axial- vector currents. They are defined to transform under local (gi, go) as 

L/,^ L'^ = g L L^ g[ + ig L d^g[ , 

R n >-* K = 9rR^ 9r + igR d^g R ■ (3.4) 

The field strengths 

Lfiv = o^Ly — o u L^ — [Lp, L u \ , R^ u = o^R u — o^R^ — [R^, R u \ , (3-5) 
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hence transform covariantly. The symmetry breaking parameters are promoted to 
the spurion fields with the appropriate transformation laws. For instance, let Ai(x) 
be the spurion field corresponding to the quark mass matrix M and generalize the 
mass term in the partially quenched QCD to the form 

-~Q R M£l L -~£l L M j £l R . (3.6) 

The spurion field is assumed to transform under the local (gi, g R ) G G as 

M^M' = g R Mg[. (3.7) 

We write up the low-energy effective Lagrangian in terms of Nambu-Goldstone boson 
fields whose generating functional of connected Green functions exhibits the same 
transformation as that of the microscopic theory ||5|, 46| . The parameters are then 



inserted at the positions compatible with the way how the chiral symmetry is broken 
by them in the Feynman diagrams in the low energy effective theory. In our context 
the ?7(l)em-charge matrix Q must also be promoted to a pair of spurion fields Ql, 
Q R which transform respectively as 

Ql^Q! l = 9lQl9 ] l , Qr^Q r = 9rQr9 r . (3.8) 

The trace of any integral multiple of Ql (Qr) is chirally invariant. It is thus possible 
to impose the chirally invariant condition 

sir (Q R ) = str (Q L ) . (3.9) 

They are not required here to vanish. After writing up the Lagrangian to the order 
of our interest, Ai is set to the diagonal quark mass matrix 

M d = diag (m\, ■■■ , m v Nv , mf, ■■■ , m s Ns , m\ , ■ • • , m v Nv ) , (3.10) 

and Ql, Q r are set to the diagonal EM charge matrix 

Q = diag (QX, ■ ■ ■ , Q v Nv , Ql , • • • , Q s Ns , QX,---, Qn v ) ■ (3.H) 

The substitution 

L M h- L M + eQA^ , R„ .-> R„ + eQA„ . (3.12) 

introduces the coupling of photons to the meson fields. 

In practice, in order to write down chirally invariant operators, it is convenient to 
use the building blocks O which transform asCn h((gL, g R )', n) O h((g L , g R ); II) 1. 
A set of the building blocks O, each of which is also the eigenstates of charge conju- 
gation and intrinsic parity transformation, is 

m m = % \u x (dyU - iR^u) - u (d^u ] - iL^u ] ) } , 

X± = u ] xu ] ±u\ ] u } x = ^B M, 

Ql = uQ l u\ Q r = u ] Q R u, 

^±^u = uL^ u u ] ±u [ R^ LU u, (3.13) 
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where Bo is the mass scale characterizing the size of the chiral condensate [ 46 1 , and 
their covariant derivatives with respect to the Maurer-Cartan form 

rV = -- {u< (d^u - iR^u) + u (d^ - iL^) } . (3.14) 

The chiral Lagrangian at the leading-order, 0(p 2 ) ~ 0(e 2 ), takes the similar 
form as in the unquenched case 

F 2 
C 2 = -f str (uX + X+) 

+e 2 Cstr (Q L Q R ) -\f, v F» v - ± {d^f . (3.15) 

Using such a normalization of the decay constant Fq that Fq ~ 90 MeV, II contains 
tc + , 7T 3 '/\/2 in its matrix element. The coefficient C parametrizes the EM correction 
induced from the short distance dynamics less than the length scale ~ l//x above 
which the present effective description is valid. Since we adopt QED^ constructed 
in Sec. |2] in case space is finite, the gauge potential A^ appearing in Eq. (|3.15|) does 
not possess the components A^t, 0). 

The (bare) intrinsic parity even Lagrangian density at the next-to-leading order 
consists of (9(p 4 )-, 0(e 2 p 2 )- and C(e 4 )-terms 

12 25 

£ 4 = J2 Lj Xj + J2 e2F o K i ©i ■ ( 3 - 16 ) 

3=0 3=1 

Here, Xj (0 = 1, • • • , 12) are the terms of 0(p 4 ) 

£ — str (u^UyU^u") , 

X 1 = (str( M X)) 2 5 

X2 = str (u^iiv) str {u^u v ) , 

£ 3 = str ((u/) 2 ) , 

X 4 = str (n M w M ) str (x+) , 

X 5 = str (u^x+) , 

^6 = (str(x+)) 2 , 

X 7 = (str(x-)) 2 , 

X 8 = - str (x+ + X-) , 

X 9 = --strfly, u v \$ + ^ v ) , 

1 
4 

in = str (n^n^ + c^ v ) 

X 12 = ~ str (xl - X-) = str ( X X J ) ■ (3.17) 



£10 = 7 str (5 r +,^S r + / - 3-, ^f) 
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L>n,i2 ar e known as Hi = L\\ H 2 = £12 i46fl . Under the condition (|3.9| ) and the 
leading-order equations of motion, 2)^ (j = 1, • • • , 25) describe all possible EM 
corrections at 0(e 2 p 2 ) or 0(e 4 ) in terms of Nambu-Goldstone boson fields. The first 
19 ones are as follows; 



2Ji 

2) 2 
2)s 
2J 4 

2J 5 

2J 6 

2> 

2J 8 

2J 9 

2Jio 
2Jn 

2Jl2 
2Jl3 
2Jl4 
2Jl5 
2Jl6 

2Ji 7 

2Jl8 
2Jl9 



^strnQ L ) 2 +(Q i j) 2 )str(// /i //H 



str (g L g fl J str («X) , 

-str [QrU^ str (Qi^J - str \Q L u^j str (q l m a 

str (Qk?v) str [QlU^ 

str 



2 / ~ \ 2 



U^W 



str {{ Q R Q L + QlQr) u.u 



^strnQ L ) 2 +(Qfl) 2 )str( V+ ) 



str^Q L Q R Jstr(x + ) , 

str |(Ql) 2 +(Qr) 2 }^ 

str((Q K Q L + Q L Q R )x+ 
strffQ K Q L -Q L Q K )x- 



istr 



V„Q R , Q fl 



u A 



V„Q L , Ql 



w 



str [VuQrVQl 



str [VuQrVQu + VuQlVQl 



e z F z str [Q R Q L 



e 2 F 2 str(Q R Q L )str( ( Q R ) + ( Q L 



e 2 F 2 ( str ([Q r )\(Ql" 



str Qru^Qru 11 + Qlu^Qlu 1 - 



str (QRUnQ L u^ 



(3.18) 
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The rest 6 ones are as follows; 

2 / „ \ 2 



2) 20 = e 2 F 2 str((g K j [Q L 
%2i = e 2 F$ sir (q r Q l Q r Q l ) , 



2) 22 = e^str [Q 



Qi 



2) 23 = e'F z i str (Q R ) str I Q fi (^Q L J I + str (^Q L J str ( Q L ( Q 
2) 2 4 = str (Qk) str (q l u^ + str (g L ) str (q r u^ 
2) 25 = str (Qr) str (QlX+I + str (q l ) str (QrX- 



ifl 



(3.19) 
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kj 
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3 
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"T ^ 


18 
19 


3 

4 
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20 


2iV 5 Z 2 - 3Z 


8 


Z 


21 


2N S Z 2 + 3Z 


9 


l 

4 


22 


-Z 2 


10 


1 , JV S 7 


23 


-8Z 2 


11 


i 

8 


24 


-Z 


12 


1 

4 


25 


-Z 
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The coefficients Lj, 



Kf- absorb the 
UV divergences that arise from the one- 
loop correction 



Lj + lj A e 
Ka + ki A e 



Kf{n) . (3.20) 



Here we employ the dimensional regular- 
ization where the spatial dimension is an- 
alytically continued to d. Denoting the 
full space-time dimension as D = d + 1 = 
4 - 2e, A e in Eq. ( |O0l ) is given by 

,2 



A r 



32tt 2 



-In 






-lE + l 

(3.21) 



Table 1: Coefficients fcj of UV divergence 
in Ki. 



The values of lj (j = 1, ■ ■ • , 12) are avail- 
able in Ref. |^7j for the generic number of 
flavors. The values of kj (j — 1, • • • , 14, 
18, 19) were computed in Ref. |33j for Ns = 



3 . Table p] lists the values of kj 2 computed using the heat kernel method and eval- 
uating the fermionic pion contribution explicitly for general Ns without setting the 
quantities in Eq. ( |3.9| ) to zero. The compuatation is performed only for Feynman 
gauge A = 1. The dimensionless quantity Z in Table |] is defined by 






2 The normalization of our kj's differs from those in Ref. [B3J such that (kj) c 



(3.22) 



-2(fej)BD for 



j = 1, • • • , 10, 12, 13, 14, 18, 19, and (k n ) 



2(fcn) 



BD- 



13 



For the sake of simplifying expressions, the superscript "R" is omitted from every 
renormalized parameter that appears in what follows. 

In the momentum space, the Feynman rule in the partially quenched chiral per- 
turbation theory in finite volume is the same as that in infinite volume. In particular, 
the form of propagators that allows us to carry out the computation efficiently is 
obtained by introducing the super-traceless component called super-?/, deriving the 
propagators in the matrix element basis $ 7 j (J, J = 1, • • • , (N S + 2N V )), and taking 
the decoupling limit of super- 7/ []38], J39| . 



The entry xu of the matrix of meson mass squared including the leading-order 
EM correction is given by 

Xi + Xj , 2e 2 C 2 , Q OQ \ 

Xu = = + —pr W ~ Qj) , (3-23) 

where xi an d qi are the eigenvalues of x\m-*m an< ^ 2 in Eq. (|3.11|) , respectively. The 
propagators are shown here in the case that the eigenvalues X(r) = Xn v +t 0- < r < 
Ng) corresponding to sea quarks as well as the mass squared Xx (x = 1, • • • , N s — 1) 
of the diagonal meson eigenstates in the sea-meson sub-sector, referred to as sea 
mesons, differ from all Xj (J = 1> ' ' ' > Ny). The propagators of bosonic mesons 
(1 < I,J } K,L < (Ny + N s ) or (1 + N v + N s ) < I, J, K, L < (2N V + N s )) are 
denoted by 

iGj. K L (p 2 ) = Jtfixe** <ll'»n*(0)) . (3.24) 

The off-diagonal meson fields have the usual form of propagators 

iG I J . K L {p 2 ) = 8 I L 8 J K ^ (I^J,K^L). (3.25) 

V ~ Xu 

The propagators of the fermionic mesons H^- = IT s v a (1 < i < Ny, 1 < j < 
(Ns + Ny)) also have simple forms 



i5* J ,* I (p a ) = y'd* a ;e*-(s* i ( a r)S t *(0) 



= <V/^ (l<i,l<Ny,l<j,k<(N s + N v )). (3.26) 

V ~ Xij 

The one-loop calculation needs the propagators of diagonal mesons only for i, j = 

f At 

P 2 ~Xx / 

(3.27) 
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Esm 
denotes the sum over all sea mesons, and 

X 

N S 

II fa* ~ ^(r)) 
r> i — r=l __ jd i 

(Xi-Xj)Y[ (Xi-Xx) 

X 

N S 



II (** - *(*)) 



K ij = 



r=l 



nam 



7V S 

n (^ ~ ^w) 



^(d) = r=\ 



nam 



x 

N S 



R 



n (^ ~ *w) 

(s) __ r =l 



nam i / 

(x* - Xx) \^i * ~ *(•) 



N S 1 1 \ 

Ex — -\ sm 1 



Xi Xx 



3.2 electromagnetic correction to meson mass in finite volume 

With the preparations done in the previous subsection, we are ready to compute the 
next-to-leading order correction m^| (L) to the off-diagonal pseudoscalar meson 
mass squared in finite volume, which can be obtained from the self-energy function 
E*-(p 2 ) within our approximation as 

m 2 ij \ A (L) = T* j ( X i j )\ L , (3-29) 

and investigate the relevance of finite size correction to the EM splitting. As recalled 



in Sec. p.2| Lorentz boost symmetry is violated in finite volume. Here mfA (L) is 



«Jl4 



defined in the rest frame, p M = (y/xij, 0). 

Before carrying on the explicit calculation further, we mention the limitation 
inherent to our approach on the ability to capture the dynamics. While the low 
energy effective theory enables us to evaluate the effect of finiteness of volume on the 
virtual quanta with low frequencies, it provides no knowledge on the effect to the 
short distance dynamics packed in the low energy constants Lj(fi) and Kj{(x). The 
finite size scaling effect on the low energy constants Lj(fi) carrying the information 
on QCD less than l//i, is insignificant so as to affect to the properties of Nambu- 
Goldstone bosons. The situation is different in QED where there is no intrinsic 
scale 1/fj, that separates long and short distances [BIJ. For this reason, though the 
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finite size scaling effect on QED is dominated by the loop contribution whose finite 
part represents the long distance physics, the low energy constants Kj(/j,) possibly 
suffer from sub-dominant but non-negligible finite size scaling effect. For the purpose 
explained in Sec. [l|, even qualitative features of the finite size scaling effect such as 
the size and sign are worthwhile to investigate. 

With the above remark in mind, we proceed to calculate m 2 -\ (L) in finite 
volume. For the purpose of (1) getting m 2 A (L) for N s = 2, 3 simultaneously, and 
(2) demonstrating that QED^ share the common UV divergent structure as QED 
in infinite volume, we describe the computation in some detail. As the probes are 
valence quarks, it suffices to put 1 < i(^)j < Np. Straightforward calculation of 
four types of Feynman diagrams yields 

M-6yi^n 2 (* 2+2 *r)GV# 2 ) 

ker d 

+ (x,-fc 2 )(Gy^ 2 )-sy^ 2 )) 



^i I A 



+ (Qij) 2 e 2 



+qijq n j G n j. J n {k 2 )} 

--Wp-d — 

2tt L d ^ r ' -k 2 

p ■ k 

+2 






-kHxij-ik+p) 2 ) 
+m 2 CtlJ , (3.30) 

represents the sum over the indices n of sea quark flavors only, T d = 

n 

f d - {0} with f d in Eq. ( pO|) , and 



QCD _ Xi + Xj 

Xij — 2 ' 

qij =q%- qj ■ (3.31) 

ttiq i ■ represents the contribution of the next-to-leading order local terms in Eqs. ( |3.17| ), 



( |3.18|) and ( P-19| ) to the pseudoscalar meson mass squared, m^; in Eq. ( [3.30 



is 



written in terms of bare Lj , Kj . The explicit form of m c i ■ will be given after renor- 
malization is performed (the sum of thq^ and vr? G i3 -| in Eq. ( |3.41|) ). We note 
that the expression in Eq. ( |3.30| ) is independent of the values of the gauge parameter 
A. 
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Eq. (|3.30| ) can be rewritten in terms of four basic functions, ( |A.3|) , (|A.20|) , ( |A.26| ) 
and (|A.30|) , defined in Appendix |A|. The function Ii(m 2 ; L) in Eq. (|A.30|) emerges in 
the QCD finite scaling effect |23[. In Appendix [A], we derive the expression written in 
terms of a Jacobi theta function for the other three functions. Applying the identities 



(2 Xi + x J )R lJ l -7;Ri=7;(x i + x 1 )R l 



1 ? (d) 3 

2 " L * 2 VAi v/ ' " '■' 



(x* + X* + Xi) V + ^V^ ^ x + ^V^ V = o Oft + *) V ' ^ 3 - 32 ) 



(Xi + 2 Xi) V - 2 ^i = 2 (» + Xi) V 

X« — Xz D a; , Xj — Xx D x " 

— ^ +— 2~ ^i =2 

to the terms in Eq. ( |3.30| ) that are written solely by I\{m 2 ; L), we get a compact 
expression 

m l\S L ) = 2A^Fj (Xj + Xj) 

x | i?./ J!^; L) + R t / h( Xj ; L) + ^Ri/hiXx, L) 

-2e 2 Z^2 {QijQinhiXin] L ) + qijq n jh(Xnj] L)} 

n 

+ ( qij ) 2 e 2 {(D- 1) I°(L) + 2 J n (Xi,; L) + 4 X ij hi{Xij\ L)} 
+m 2 Cjij . (3.33) 

From Eq. (|3.33|) , the finite size scaling correction is obtained as 

AmJ(L) = mH W - m ll (°°) > ( 3 ' 34 ) 

using the quantity w^L (oo) evaluated directly in infinite volume. The use of the 
results in Appendix [A] immediately leads 

AmJ(L) = Am 2 QCD]ii (L) + Am 2 EMtij (L) , (3.35) 

where Amk CB ^(L) and Am| M ^(L) represent the finite size scaling correction on 
QCD and QED, respectively 

Am 2 QCDiii (L) = * ?|±j| -L (VM#) + Ri/Miy/xjL) 



Esm 



2e 2 Z 1 



Am EM,ij( L "> = ~7^2j^Yl (fti linM(^X^L) + q tj q nj M(^UjL) 



-3 



{qijfe 2 K 



4vr L 2 
^{M-4^^}. ,3,6, 



-17- 



This expression is free from UV divergence, confirming the statement at the end of 
Sec. [2.2[ Carrying out the renormalization for the part in infinite volume with help 
of the values of kj's in Table [I], we get the finite expression of the next-to-leading 
correction m% \ AL) 



m l\A L "> = m QCD,y| 4 ( L ) + m EM,*;l 4 ( L )' ( 3 - 37 ) 



with m 2 QCDlj \ 4 (L) (m| M>ii | 4 (L)) composed of the corresponding quantity m 2 QCDij \ A (oo) 
( m EM,ij| 4 (°°)) i n infinite volume and the finite size correction in Eq. ( |3.36| ) 



m QCD,i;| 4 (L) = m QCDf y| 4 (oo) + Am QCDi] (L) , 

l EM,y| 4 (-^) = m EM,ij| 4 (°°) + ^ m EM,Jj 



™EM,»L ( L ) = m lu,ijl (OO) + Am EM,»'( L ) • (3-38) 



Each of m QCD,iiL (°°) an d m lM,«j| 4 (°°) consists of two parts 



2 I / \ 2 / \ 1 1°°P i 2 I QCD 

m QCD, 4 i| 4 (°°) = m QCD,ii(°°)| + m C,ii| 

2 I / \ 2 / \ I loop 2 I EM /o on\ 

m EM,ii| 4 (°°) = m EM,y(°°)| + m C,ij| ' (3-39) 



The terms ^QCD,ii(°°)| > m EM,ij(°°)\ are those involving the chiral logarithms 



m QCD,«Mi loop = xtfwM r' /x ' ln (P) +R " ix > hi {js 

2 / Nlloop 2e 2 Z^S f /Xm\ . . /XnA 

m EM, y (OO) I = - ——2 2^ I «tf*n Xtn In I -J I + QijQnj Xnj In I —j- 1 

-W^Mf)- 4 }' (3 - 40) 



while the terms m^ $,- ' and m^ i A are those given in terms of low energy 
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constants 



m 



iQCD 



C,ij\ 



I EM 



m 



C,ij\ 



-$ [-4 {2 L 4 N s xs + L 5 (Xi + Xj)} Xij 

+8 L 6 N s x s (Xi + Xj) + 4 L 8 ( Xi + Xi) 2 ] , 
{-4iW (id + K 2 ) - 4 (g* + gj) (# 5 + K 6 ) 

-4 g^- (2^i8 + #19) - 4iV s Q ( ft + qj) K 24 } Xij 
+2NsW(Xi + Xj)K 7 

+2N S [W(Xi + X j) + 2x 5 (ft ~ Qjf) Ks 

+4 (q-Xi + Q 2 Xj) K 9 

+4 {qhi + QjXj + (Qi -Qjf (Xi + Xj)} K 10 

-4 (qt - qj) 2 (Xi + Xj) K n + AN S Q (q t Xi + QjXj) K 2z] 



+e% 2 



4NsQ 2 



qj ) 2 (K 15 + K w ) + 2(q*-q*) 2 K: 
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+4( gi - q 3 f (q 2 + gj) K 21 
+4N s Q(qi -qjf (q t + qj) K 23 ] . 



(3.41) 



In the above, 



N S 



1 N S 

at: z^ X{r) ' 

15 r=l 



_ 1 ^S ^ N S 

a r=l ° r=l 



(3.42) 



<7iV- 



v+r- 



The formulas derived thus far 



with use of electric sea quark charges g( r ) 

applies both to Ns = 3 and Ns = 2 if the sum over sea flavors and the one over 



sea mesons are appropriately understood. The expressions for mQ CD 

2 



*?U 



m 



(oo) and 
EM,i?| 4 (°°) wm re duce to the result found in Ref. |33[ if one sets N s = 3 and 
Q = and discards the terms of order e 4 . 

Before turning to the numerical analysis, we observe a few features that can be 
read off from Eq. (|3.36|) . First, the asymptotic behavior for L^xJj ^ 1 is determined 
by the term including the function TC(x). This term behaves like 1/L times the 
function Ti^XijL)- As can be seen from Fig. p], TC(x) gradually increases for x —>■ oo. 
We find that the decrease of the finite size scaling effect is slightly slower than 1/L. 

Secondly, the terms in Eq. (|3.36|) all vanish for y/xu ~^ except for the one 
proportional to the constant k defined in Eq. (|A.29|) 3 . This term seems to appear 
whatever is used as the low effective field theories, because the contribution lead- 
ing to this term originates from the diagrams in scalar QED theory. For instance, 



3 Hcre we use the fact that 



H{x) 



is finite for all x. This point, however, has been checked only 



through our numerical investigation which gives Eq. (A. 19) 
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that contribution is also involved in the model including the vector and axial- vector 
mesons that leads Eq. ( |1.2|) , as is suggested from the overall numerical factor. The 
presence of this term indicates that 1/L should be regarded as being the same order 
as the pseudo-Goldstone mass and the elementary charge e in order for the chiral 
perturbation to remain systematic. In fact, the calculation done thus far implicitly 
assumes the relative magnitude of pseudo-Goldstone mass and L in the p-regime 
[Pm, m n ~ ~ ~ p. Our result suggests that there is a p-regime for the systematics 
chiral perturbation theory including electromagnetism in finite volume. The issue 
examining whether this is actually true is beyond the scope of this paper. 

3.3 numerical investigation 

We turn to the numerical evaluation of the EM correction in the next-to-leading 
order approximation 



2e 2 C 

-Tfl 



2 / r \ Ze ° / \2 . 2 



"'i^r .-,■(£) = —EvT (ft ~ 1i) + m EM.ij\ 4 ( L ) , (3-43) 





to study the qualitative features of the QED finite size scaling effect. For that purpose 



the values of various low energy constants found in Ref. |58| are used as reference 



F = 87.7 MeV, C = 4.2 • 1(T S GeV 4 , 

L 4 = 0, L 5 = 0.97 ■ 10~ 3 , L 6 = 0, L 8 = 0.60 • 10~ 3 , 

K 5 = 2.85 ■ 10~ 3 , K 9 = 1.3 ■ 10" 3 , K w = 4.0 ■ lO" 3 , 

K n = -1.25-10" 3 . (3.44) 

The others are set to zero. We employ the formula for iVg = 3 and Ny = 3 and set 

Xs = X6 = X9 = (500 MeV) 2 , 

2 1 

q x = g 4 = q 7 = -, g 2 = q 3 = q 5 = q 6 = q 8 = q g = -- , (3.45) 

throughout the analysis. In what follows, it is always understood that the value of 
the mass of a ghost quark is set equal to that of the valence quark with the same 
flavor. 

We plot the dependence of m| M 12 (L)/mg M 12 (oo) on L in Fig. |2| for the val- 
ues of quark masses corresponding to Xi = X2 = Xa = Xz = (150 MeV) 2 and 
(300 MeV) 2 . The horizontal axis denotes the linear size of volume L normalized in 
unit of 1/Mp ~ 1/(770 MeV). For instance, for the size L = 16a ~ 0.72 x 1/M p 
used in Ref. @, m 2 M 12 (L)/m 2 M12 (oo) are 0.144 and 0.421 for X i = (150 MeV) 2 
and Xi = (300 MeV) 2 , respectively. Thus, our calculation indicates that the finite 
size effect is significant for an available lattice geometry. The difference between two 
quark masses emerges for small L, and the finite size effect is smaller for larger quark 
mass. 
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We next study the EM splitting Amf^L) in the (valence) kaon mass squared. 
Using m| M :ij (L) in Eq. ( [3.431 ), it is given by 

Am|(L) = m| M13 (L) - m| M23 (L) . (3.46) 

Figure |3| shows the L-dependence of Am|^(L)/Am^(oo) for the two sets of quark 
masses corresponding to the same values of xi — X2 = Xi — Xs used in Fig. |] with 
X3 = (500 MeV) 2 fixed. The L-dependence in Fig. ^ is almost similar to that in Fig. |2] 
for each set of masses. For small L and Xi — (150 MeV) 2 , the size of finite size scaling 
effect is smaller than the electromagnetic correction to the charged pion mass. For 
instance, for L = 16 a, Am|(L) / Am|(oo) is 0.299 for xi = (150 MeV) 2 and 0.415 for 
Xi = (500 MeV) 2 . Figures |] and || show that the values of EM correction depends on 
the quark masses for L = 16 a. This observation indicates that the sizes of the terms 
depending on the quark masses in Eq. ( |3.36| ) are comparable to or more important 
than that of the term proportional to k for the quark masses used in the present 
analysis. 

Figure |] compares the L-dependence of the electromagnetic splitting of m 2 K (L) 
in partially quenched QCD with that in (unquenched) QCD. The two sets of plots 
are drawn in Fig. [| for 

set 1 <S> Xi = X2 = (150 MeV) 2 , X a = Xs = (300 MeV) 2 , 
set 2 <£> xi = X2 = (300 MeV) 2 , X 4 = Xs = (150 MeV) 2 . (3.47) 

As can be seen, these two sets of plots in partially quenched QCD coincide with each 
other. By evaluating the L-dependence for various set of the values of quark masses, 
we find that the L-dependence for a set of the values of quark masses is almost the 
same as that for the set with the valence and sea quark masses interchanged. This 
observation together with the comparison of the relative sizes of four sets of plots 
in Fig. |] shows that the size of the finite size correction is roughly determined by 
the average value of the quark masses involved irrespective of whether the system is 
unquenched or partially-quenched. To elucidate this point, we show in Figs. |5] and |6] 
the quark mass dependence of the finite size scaling effect on m 2 K (L) on m? = xi — 
X2 = ^BQin u for fixed L = 16a and L = 32a (a ~ 1/(1.66 GeV) [[L2"|), respectively. 
In the partially quenched case, the masses of up and down sea quarks are fixed to 
be 300 MeV. We can see that, at L = 16 a, the size of the finite size effect changes 
rapidly for m 2 < 0.05 [GeV] 2 in the unquenched case, while no such significant 
change is observed in the partially quenched case. For larger m 2 , they approach to 
each other. For L = 32 a and larger L, the finite size effect is determined by the 
relative size of L and X3 — Xq — (500 MeV) 2 . 

4. Conclusion and discussion 

In this paper, we studied the finite size scaling effect on the electromagnetic (EM) 
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correction to the pseudoscalar meson masses from the low energy effective field theory 
of QCD including electromagnetism. For that purpose, we began with constructing 
a new QED in finite volume, QED^. Taking the practical application to the lattice 
simulation into account, we adapted it to the partially quenched chiral perturbation 
theory including electromagnetism. We computed the electromagnetic correction 
to the pseudoscalar meson mass squared at the next-to-leading order on both of the 
spaces with finite and infinite volumes for generic number N$ of sea quarks. Through 
numerical investigation for N$ = 3, we found that the finite size scaling effect on 
the EM correction is sizable on the space with the volume available in the lattice 
simulation. By investigating its dependence on the quark masses in unquenched and 
partially quenched systems, we pointed out that the finite size correction is deter- 
mined by the averaged values of masses of quarks involved in the system. Though 
the current study was restricted to the pseudoscalar meson masses, we can study the 
EM corrections to the other hadronic observables such as decay constants in finite 
volume. 

It should be noted that the finite size correction tends to increase for the increase 
of the extent T of the temporal direction. For example, let us suppose that the value 
~ 0.9 in Eq. ( |1.4|) does not change so drastically even if the contribution of A^(t, 0) is 
subtracted. The explicit calculation in the low energy effective theory including vec- 
tor and axial- vector mesons show that it is reduced to about 0.6 in the limit T — > oo. 
Since the simulation has to be performed with finite T, m| M 12 (L)/mg M 12 (oo) ob- 
tained in our present analysis with T — > oo may overestimate the finite size correction 
in the actual simulation. It is one of the subjects to fill the gap between these two. 

Finally, we recall that we construct QED^ because we would like to respect 
the boundary condition (i.e., periodic boundary condition) for the meson fields along 
spatial directions usually assumed in the lattice QCD simulations, and to understand 
the qualitative properties of the finite size scaling effect with the practical situation 
in our mind. It is plausible that the finite size scaling effect depends on the boundary 
conditions imposed on fields. As illustrated in Figs. f|, if periodic boundary condition 
is employed, the dependence on the volume is not simple around the available volume 
L ■ (0.77 GeV) < 10 in the lattice simulations. It would be one of the important 
subjects to study the finite size scaling effect on EM corrections for the various types 
of boundary conditions and to select the one in which the finite size effect exhibits 
such a simplest behavior that allows us to extrapolate the data at available volumes 
to the value in infinite volume. 
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A. Formulae for several sums 

The one-loop corrections to meson masses in chiral perturbation theory for QCD 
plus QED system with finite volume are described in terms of several basic functions. 
Each of these functions takes the form of the one-dimensional integral of the sum over 
three-dimensional momenta of some function. We employ dimensional regularization 
and define 

°° dk° 1 v-^ , A . 

ker d 
Here, d = D — 1, V = L d and the sum runs over all k e Yd, where 

f d = jk=(fc\ ■■■ ,k d ) W e^zX . (A.2) 

The aim of this Appendix is to get compact expressions for the four functions that can 
be evaluated by MATHEMATICA and so forth by following the strategy in Ref. [p3" ]. 
We first consider the function 



, , 2 „ / i-a-i f° dk l v- 



1 



^ 2m V ^ (-P - ie) {m 2 - (k + p) 2 - ie) 



. (A.3) 



In the above, p is assumed to be on-shell, p 2 = m 2 and T' d = Td — {0}. 
By introducing the Feynman parameters as usual, Eq. (|A.3j) becomes 



k 6 r d 



1 



2m V ^ {y 2 m 2 - (k + yp) 2 - ie} 2 ' 
ker d 

The sum of a function F(1&) over k G T' d can be written as an integral 

v E ^ (k) = / ddk '^) \ E 6d ( k ' - k ) - v ^ (0) • 



(A.4) 



V 

ker' ker d 



(A.5) 



Using the Poisson resummation 



7E^-') = ^E^ < A - 6 > 



ker d v ; xer 

where 



r d = {x= (x 1 , ••• ,x d ) \x j eLZ} , (A.7) 
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the sum (IA.5I) can be written as 



ker' x e r d 



(A. 



where -F(x) is obtained by Fourier transformation of -F(k) 

d d k 



F(x) 



ikx 



(2tt)< 



F(k) 



(A.9) 



In the present context -F(k) is 



F(k) 



1 



{y 2 m 2 — (k + yp) — ie} 

rM / T ^^ exp ^~ iA ^ m2 _ ^ + 2 ^ 2 ~ ie ^ 



(A.10) 



A given (i-dimensional vector x defines a .D-dimensional space-like vector x^ = 
(0, x). A one-dimensional integral in Eq. ( |A.4| ) and a d-dimension integral in Eq. (|A.9| ) 
combine to become an integral over .D-dimensional momenta, which can be readily 
carried out; 



~r ( iA ) 2 / ~tttd ex P [~ iX {-( k + yp) 2 ~ ie } 

A J l (Z7r) 



IX 



*] 



dX X 2 



— exp 



~ A 1 v^x ~ yp ) ~ ( yp ^ 



'0 A (4ttA) 
For p^ = (m, 0), [i-^ — yp) — (yp) 2 = jg*. Eq. ( |A.4| ) thus becomes 



(A.11; 



I u (m 2 ; L) = (^ 2 ) 



2\ 2 -- 



jf*(S-/ A ) 



X 



00 dX 

T 



(4ttA) 



D" eX P 

2 



7~ 2 

\ 2 2 "^ I |2 

-Xy m n 

y 4A ' ' 



In the sum appearing above, the contribution of n = is exactly Iu(m 2 ) 
00). Iu(m 2 ] L) — Iuirn 2 ) is hence free of UV divergence. Letting D 
difference and rescaling A — > j^X leads 



/ii(m 2 ;L)-/n(m 2 ) 



1 1 



167T 2 rnL 



H{mL) . 



(A.12) 

Iuirn 2 ; L - 
> 4 for this 



(A.13) 
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Here 



H{mL) =n J ^ erf ( mii/^- ) S( A ) ( A. I I ) 



S(X) = -[ Y, -/^]exp(-^|n| 2 ) 

,ne23-{0} J J 



1 \ \ 3 



t? 3 (0, ij) ] -1-A^, (A.15) 



where erf (x) is the error function 



erf(x) = -%= f dse~ s \ (A.16) 



o 



and $3(u; r) is a Jacobi-theta function 

oo 

173(1;; r) = 2_, exp (^rm 2 + 27Tf m) . (A. 17) 

n=— oo 

We recall that the term A2 in 5(A) arises in our new QED, QED L . Because 

^ 3 (°'4) ^ a1, (A - 18) 

in the infrared limit A — > 00, the presence of that term indeed ensures IR-fmiteness 
of the integral over A in the expression (|A.14 ) . The numerical study shows that 



l im ZlM ~ 10.4. (A.19) 

x-+0 X 



Next we consider the function 

2 v 2 -f f°° dk° 1 ^ p-k 



T / 2 r x / 2^2-^ f°° dk" 1 v^ 

j u (m»iL) = o?) ■jf_ SS7 E 

(A.20) 



27Ti V ^ (-P _ i£ ) | m 2 _ ( fc + p )2 _ ^J 

ker d 



The same manipulation as done for I n (m , L) in Eq. ( |A.10| ) leads 



Jn(m 2 ; 



^-r*^) 



X 



1 r dx c\\*( *\*-% f 



r(2)J A ^"' v~/ y z(27r) D 

p ■ kexp [— (iA) |i/ 2 m 2 — (/e + i/p) — ie} — ix • fe] . (A. 21) 
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We follow Appendix of Ref. [4£] to carry out such integrals with k^ in the numerator. 
We replace k^ appearing in the numerator as 



Ki, 



1 1 d 



exp (—2 (iX) p ■ k) 



Jp^O 



2 iX dpi 1 

Carrying out the integral of k and performing a p-derivative gives 

00 dX X 2 



(A.22) 



Jn (m ; 



l) ^H?^/ a H 



XP ' * l 2A ~ VP ' 6XP 



£(£4/* 



2 \ ^ V 



.r 



A (4ttA) t 

D dA A 

o ^ (4ttA)"^ 



W>\ l x l 2 

1 - e" Am J e « , 



(A.23) 



where the second equality follows by taking p^ = (m, 0) and performing the integral 
over y. In Eq. ( |A.23| ), UV divergence is contained in the term with x = 0, which is 
exactly Jn(m 2 ) = J u (m 2 , L — > oo). Therefore, Jn(m 2 ; L) — Ju(m 2 ) is UV-finite. 
In terms of the function 



it can be expressed as 

Jn(m 2 ; L) - Jn(m 2 ) 

The quantity 

dk° 1 



1 1 

32^2 J2 



JC(mL) . 



I°i(L) ee (/x 2 )' 



2?riy 



Z^ _fc2 



ker' 



—A; 2 — ie 



(A.24) 



(A.25) 



(A.26) 



needs a special care. In dimensional regularization, we put 

d D k 1 



J 1 °(oo) = ( / , 2 ) 2 ^| 



z(2tt 



,-D 



-A; 2 — ie 



However, we cannot set Ii(L) = 0. For instance, if the integral 



2\2- 



(^ 2 ) 



dk° 1 

27^17 



y _J 



ker' 



nr 



w-*r^E 



fc 2 — ie) (m 2 — k 2 — ie) 
1 



27rz L d -^ V —A; 2 — ie m 2 — k 2 — ie 
ker' 



(A.27) 
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were evaluated in two different ways; (1) direct evaluation of the left-hand side by 
introducing a Feynman parameter, and (2) evaluation of the right-hand only with 
the second term kept, inconsistency would arise. A straightforward calculation yields 



i?(L)-i?(oo) 



K 1 

AnT 2 



where k is a constant defined by 

k= I ^S(A)^ 2.837. 



(A.28) 



(A.29) 



All the one-loop contributions induced by quartic couplings are described by a 
function 



h(m 2 ; L) = (/i 2 ) : 



00 ,vz~o 



dk° 1 



£ 



2iri V *—? m 2 — k 2 — ie 
ker d 



(A.30) 



The expression for Ii(m 2 ; L) in terms of Jacobi-theta function was obtained in 

1 M{mL) 



Ref. 23 



him 2 ; L)-h(m 2 



(4vr) 2 L 2 



where 



and 



Mix) = An I ^exp(-^A)T(A), 



T ( A ) - Yl exp (~~ ' n ' 2 



(AM) 



(A.32) 



nG(Z 3 -{0}) 



(-)■,, lQ,ij 



A 



1. 



(A.33) 
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Figure 1: profile of function TL{x) 
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Figure 2: Linear volume size (L) dependence of the electromagnetic correction in the 
charged pion mass squared 
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Figure 3: Linear volume size (L) dependence of the electromagnetic splitting (Am^-) in 
the kaon mass squared 
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Figure 4: Finite size scaling effect on the EM splitting (Am^-) in kaon mass squared in 
partially quenched QCD, where open squares stand for the L-dependence for \i = Xi = 
150 MeV, X4 = X5 = 300 MeV, and dark squares with cross marks stand for the one for 
Xi = X2 = 300 MeV, X4 = X5 = 150 MeV. No significant change is observed between 
them. They differ from the other plots corresponding to unquenched QCD in the small 
volume region. 
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Figure 5: Dependence of the finite size scaling effect on the electromagnetic splitting 
(Am|) in kaon mass squared on m 2 = xi at the tree level for L = 16 a. The dark circular 
dots represent the unquenched case while the inverted triangles represent the valence quark 
mass dependence in partially quenched QCD. 
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